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This article is an overview of recent experimental and theoretical work on transport in phase- 
coherent hybrid nanostructures, with particular emphasis on dc electrical conduction. A sum- 
mary of multiple scattering theory and the quasi-classical methods is presented and comparisons 
between the two are made. Several paradigms of phase-coherent transport are discussed, includ- 
ing zero-bias anomalies, reentrant and long range proximity effects, Andreev interferometers and 
superconductivity-induced conductance suppression. 
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I. INTRODUCTION 



During the past few years, phase-coherent quasi-particlc transport in hybrid superconducting structures has emerged 
as a new field of study, bringing together the hitherto separate areas of superconductivity and mesoscopic physics 
[Altshuler, Lee and Webb 1991, Buot 1993]. The seeds of this new field of research were sown for a variety of different 
reasons. The search for novel devices led several groups to embark on experimental programmes aimed at exploring 
the properties of hybrid semiconducting-superconducting structures. During the 1980s there had been great progress 
in understanding transport properties of normal sub-micron conductors and therefore towards the end of that decade, 
it seemed natural to ask how these were affected by the presence of superconductivity. Partly as a consequence of 
formidable technological problems to be overcome in growing hybrid nanostructures, the first experimental results 
did not emerge until the summer of 1991, when two groups [Petrashov and Antonov 1991, Kastalskii et al 1991] 
independently submitted papers on phase-coherent transport in sub-micron superconducting structures. The first of 
these reported an increase in resistance in metallic conductors due to the onset of superconductivity and the second 
the discovery of a zero-bias anomaly in the dc resistance of a superconductor - semiconductor (S-Sm) contact. The 
former is still the subject of theoretical investigation [Wilhelm et al 1997, Seviour et al 1997], while the latter is now 
well- understood. 

In a phase-coherent normal-superconducting (N-S) structure, the phase of quasi-particles as well as Cooper pairs is 
preserved and transport properties depend in detail on the quasi-particle states produced by elastic scattering from 
inhomogeneities and boundaries. A key physical phenomenon, which arises in the presence of superconductivity is the 
possibility that an electron can coherently evolve into a hole and vice versa. This phenomenon, known as Andreev 
scattering [Andreev 1964], occurs without phase breaking and is describable by a variety of theoretical techniques. The 
effect of superconductivity on transport across a N-S interface is of course an old subject. In lowest order, the classical 
tunneling Hamiltonian approach ignores Andreev scattering and predicts that the dc conductance G is proportional 
to the density of states. Later [Shclankov 1980, Blonder, Tinkham and Klapwijk (BTK) 1982, Blonder and Tinkham 
1983, Shclankov 1984] it was pointed out that the contribution to the sub-gap conductance from Andreev scattering 
can be significant and a theory of a clean N-I-S interface was developed, which showed that for a delta-function 
barrier, there is indeed a marked deviation from tunneling theory, but as the barrier strength is increased, the result 
of classical tunneling theory is recovered. BTK theory applies to a one-dimensional N-I-S system or, by summing over 
all transverse wavevectors, to 2 or 3 dimensional systems with translational invariance in the plane of the barrier and 
yields for the current I through the contact 

/oo 
dE (f(E - eV) - j{E)){l + A{E) - B(E)) (1) 
-oo 

where A(E) and B(E) are Andreev and normal reflection coefficients listed in table II of [Blonder et al 1982], f(E) 
the Fermi function and fl a measure of the area of the junction. In the presence of disorder or other inhomogeneities, 
this must be replaced by the more general expressions outlined in sections II and IV below. 

Prior to 1991, experiments on N-I-S point contacts have been in broad agreement with BTK theory, exhibiting a 
conductance minimum at zero voltage V = and a peak at eV ~ A, where A is the superconducting energy gap. 
In the experiment of Kastalskii et al. [1991], the dc current through a Nb-InGaAs contact is measured as function 
of the applied voltage. At the interface, depending on the semiconductor doping level, a Schottky barrier naturally 
forms so that the system behaves like a superconductor-insulator-normal (S-I-N) structure. (An exception to this is 
InAs, which does not form a Schottky barrier at an N-S interface.) According to BTK theory, as the barrier strength 
increases the sub-gap conductance should vanish. In contrast, the experiment revealed an excess sub-gap conductance 
at low bias, whose value was comparable with the conductance arising when the superconducting electrode is in the 
normal state. This zero bias anomaly (ZBA) was later observed by Nguyen et al. [1992] in an experiment involving 
InAs-AlSb quantum wells attached to superconducting Nb contacts and by using high transmittance Nb-Ag (or Al) 
contacts of varying geometry, Xiong et al. [1993] were able to observe the evolution from BTK to ZBA behaviour. 
In an experiment by Bakker et al. [1994] involving a silicon-based two-dimensional electron gas (2DEG) contacted to 
two superconducting electrodes, a gate voltage was also used to control the strength of the ZBA and in [Magnee et al 
1994] an extensive study of the ZBA in Nb/Si structures was performed. Since these early experiments, a great deal 
of effort has been aimed at observing Andreev scattering in ballistic 2DEGs, including [van Wees et al 1994, Dimoulas 
et al 1995, Marsh et al 1994, Takayanagi and Akazaki 1995(a), Takayanagi, Toyoda and Akazaki 1996(a)] 

Kastalskii et al [1991] attributed the excess conductance to a non-equilibrium proximity effect, in which supercon- 
ductivity is induced in the normal electrode, giving rise to an excess pair current. Initially this phenomenon was 
seen as separate from Andreev reflection, but subsequent theoretical developments have shown that the distinction 
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between the proximity effect and Andreev reflection is artificial. As will become clear later, the ZBA arises from an 
interplay between Andreev scattering and disorder-induced scattering in the normal electrode. Andreev scattering is 
sensitive to the breaking of time reversal symmetry and as a consequence the conductance peak is destroyed by the 
introduction of a magnetic field. 

Zero bias anomalies constitute the first of a small number of paradigms of phase-coherent transport in hybrid N-S 
structures. A second paradigm is the observation of re-entrant [van Wees et al 1994, Charlat et al 1996] and long-range 
[Courtois et al 1996 ] behaviour signalled by the appearance of finite-bias anomalies (FBAs) in the conductance of 
high-quality N-S interfaces. At high temperatures T > T* and bias- voltages V > V* , where for a N- metal of length L 
and diffusion coefficient D, ksT* = eV* = \f(D/L 2 ), both the ZBA and FBA conductance peaks decay as l/VT and 
\j\/V . For a clean interface there also exists a conductance maximum at V*, T* and therefore at low-temperature and 
voltage a re-entrance to the low-conductance state occurs. An interesting feature of this phenomenon is the long-range 
nature of the effect, which typically decays as a power-law in L*/L, where L* = y^D/eV). This behaviour is in sharp 
contrast with the exponential decay of the Josephson effect and has been observed in a number of experiments. In a 
T-shaped Ag sample with Al islands at different distances from the current-voltage probes [Petrashov et al 1993(b)], 
a long-range proximity effect was observed, in which the influence of the island extended over length scales greater 
than the thermal coherence length L*. Similar behaviour was also observed [Petrashov et al 1994] in ferromagetic- 
superconductor hybrids made from Ni-Sn and Ni-Pb. In an experiment involving a square Cu loop in contact with 2 
Al electrodes Courtois et al [1996] clearly identified both the short- and long-range contributions to phase-coherent 
transport. In this interferometer experiment, they observed a phase-periodic conductance decaying as a power-law in 
1/T, in parallel with a Josephson current which decays exponentially with L/Lt- 

The above re-entrance phenomenon is also observed in a third paradigm of phase-coherent transport, which arises 
when a normal metal is in contact with two superconductors, with order parameters phases 4>i and cf>2, whose difference 
cf> = (pi — <f>2 can be varied by some external means. Prior to the experimental realisation of these structures, the 
electrical conductance of such Andreev interferometers was predicted to be an oscillatory function of cf). Spivak and 
Khmcl'mtskii [1982] and Al'tshuler and Spivak [1987] identified a high temperature (T >> T*), weak localisation 
contribution to the conductance of a disordered sample, whose amplitude of oscillation was less than or of order 2e 2 /h. 
For an individual sample, the period of oscillation was found to be 2tt, but for the ensemble average a period of 7r 
was predicted. Nakano and Takayanagi [1991] and Takagi [1992] examined a clean interferometer in one-dimension 
and again predicted a 27r-periodic conductance with an amplitude of oscillation less than or of order 2e 2 /h. Lambert 
[1993] examined a disordered conductor in the low-temperature limit (T << T*) and identified a new contribution 
to the ensemble averaged conductance with a periodicity of 2ir. This 2ir periodicity is a consequence of particle-hole 
symmetry, which also guarantees that at zero temperature and voltage, the conductance should possess a zero phase 
cxtrcmum [Lambert 1994]. Prior to experiments on such devices, the generic nature of this prediction was confirmed 
in numerical simulations [Hui and Lambert 1993(a)] encompassing the ballistic, diffusive and almost localised regimes 
and in a tunnelling calculation of the ensemble averaged conductance by Hekking and Nazarov [1993]. 

The first experimental realisations of Andreev interferometers came almost simultaneously from three separate 
groups. In March of 1994, de Vegvar et al [1994] showed results for a structure formed from two Nb electrodes in 
contact with an Al wire. They found a small oscillation 10~ 3 (2e 2 //i) with a sample specific phase in the 27r periodic 
component, suggesting that the ensemble averaged conductance should have a periodicity of 7r, in agreement with 
Spivak and Khmcl'nitskii [1982]. However in contrast with all subsequent experiments, no zero phase extremum was 
observed. In April/May of that year, Pothier et al [1994] produced an interferometer involving two tunnel junctions, 
which showed a 27r-periodic conductance, with a zero phase maximum and a low-bias, low-temperature amplitude of 
oscillation of order 10~ 2 (2e 2 //i), which decayed with increasing temperature. In May 1994, van Wees et al [1994], [see 
also Dimoulas et al 1995] produced the first quasi-ballistic InAs 2DEG interferometer, with high transparency N-S 
interfaces. This experiment showed the first re-entrant behaviour in which the amplitude of oscillation SG varied from 
SG « — 0.08(2e 2 //i) at zero voltage, ( where a minus sign indicates a zero phase minimum and a + sign a zero-phase 
maximum) passes through zero at a bias of order O.lmV, reaches a maximum at a bias of order V* and then decays 
to zero at higher voltages. (For a detailed study see [den Hartog et al, 1996]). Unlike the Josephson current which 
decays exponentially with T/T*, these conductance oscillations decayed only as a power- law. 

The first experiment showing an amplitude of oscillation greater than 2e 2 /h was carried out by Petrashov et al 
[1995]. Here, silver or antimony wires in the shape of a cross, make two separate contacts with superconducting Al 
and the phase difference between the contacts is varied using either an external field applied to a superconducting 
loop or by passing a supercurrent throught the Al. The amplitude was found to be 5G « 100(2e 2 /ft.) for Ag and 
3.10~ 2 (2e 2 //i) for Sb, and exhibited a periodicity of 2n. In this experiment, the phase difference <j) was varied both by 
passing a magnetic flux through an external superconducitng loop and by passing a supercurrent through a straight 
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section of the superconductor, thereby emphasising that precise the manner in which the order parameter phase is 
controlled is not important. These experiments were crucial in demonstrating that in metallic samples, the ensemble 
averaged conductance is the relevant quantity and therefore a quasi-classical description is relevant. It is perhaps 
worth mentioning that with hindsight, an earlier experiment reporting large-scale oscillations in a sample with two 
superconducting islands [Petrashov et al 1993(a)] can be regarded a precursor to these interferometer experiments. 
However the phase of the islands was not explicitly controlled, making an interpretation of the 1993 experiments more 
difficult. 

A fourth and more recent paradigm of phase-coherent transport is the appearance of negative multi-probe conduc- 
tances in structures where Andreev transmission of quasi-particles is a dominant process [Allsopp et al 1994]. The 
first experiments reporting this behaviour were carried out by Hartog et al [1996], using a diffusive InAs 2DEG. These 
probe individual coefficients in the current- voltage relations and demonstrate fundamental reciprocity relations arising 
from time-reversal and particle-hole symmetry. 

Finally a fifth paradigm is the suppression of electrical conductance by superconductivity in metallic systems 
without tunnel barriers. Experimentally this phenomenon has been observed in several structures, involving both 
non-magnetic (ie Silver) and magnetic (ie Nickel) N-components in contact a superconductor via clean interfaces. 
In the experiment of Petrashov and Antonov [1991] the conductance of a Ag wire on which several Pb islands are 
deposited, decreases by 3% when the islands become superconducting, this corresponding to a conductance decrease 
of SG w 100 (2e 2 /h). These experiments reveal two regimes: as the field decreases below Hc 2 , an initial resistance 
increase occurs. The resistance remains field insensistive until a lower field is reached (corresponding to the order of 
a flux quantum through the sample), at which point a further conductance suppression occurs. In the experiment by 
Petrashov et al [1993b] involving a T-shaped Ag sample with Al islands at different distances from the current-voltage 
probes, a superconductivity induced change in the resistance by up to 30% was observed, but it was found that for 
different samples, the resistance could either increase or decrease. Although superconductivity-induced conductance 
suppression in metallic samples was predicted a number of years ago [Hui and Lambert, 1993(b)] these experiments 
at first sight appeared to conflict with quasi-classical theories, which universally predict that the normal-state, zero- 
temperature, zero-bias conductance Gat is identical the conductance Gns in the superconducting state. This effect 
is addressed in [Hui and Lambert 1993b, Claughton et al 1995, Wilhelm et al 1997, Seviour et al 1997]. 

The main aim of this review is to outline the quasi-classical and multiple scattering theories needed to describe the 
dc electrical conductance of phase-coherent hybrid N-S structures. For this reason we shall not discuss thermodynamic 
phenomena such as the Josephson effect in any detail, despite the fact that ground-breaking experiments using clean 
[Takayanagi and Akazaki, 1995(b,c,d)] superconducting quantum S-2DEG-S point contacts show quantization of the 
critical current as predicted by Furasaki et al [1991] and for shorter junctions by Bcenaker and van Houtcn [1991]. 
There are several notable theoretical papers addressing Andreev scattering in such structures, including [van Wees 
et al 1991, Bagwell 1992, Furasaki et al 1994, Gusenheimer and Zaikin 1994, Zyuzin 1994, Hurd and Wendin 1994 
and 1995, Bratus et al 1995, Chang et al 1995, Koyama, Takane and Ebisawa 1995 and 1996, Levy Yeyati et al 1996, 
Martin- Rodero 1996, Wendin and Shumeiko 1996, Reidel et al 1996, Volkov and Takayanagi 1996(b)]. Similarly to 
restrict the length of this review, recent theories of thermoelectric coefficients [ Bagwell and Alam 1992, Claughton 
and Lambert 1996] and shot noise will not be discussed [Datta 1995], nor will we discuss work on coulomb effects in 
superconducting islands [Eilcs et al 1993, Lafarge et al 1993, Tuominen et al 1992 and 1993, Hergenrother et al 1994, 
Black et al 1996, Hckking et al 1993]. 

While Andreev interference effects are generic phenomena, their manifestation in a given experiment is sensitive to 
many parameters. For the purpose of writing this review, it is therefore convenient to adopt a simple classification of 
experimental arrangements sketched in figure 1. The generic structure shown in figure la represents our first class of 
N-S-N hybrids and has many realisations. For convenience, we distinguish these from a second class of N-S hybrids 
of the kind shown in figure lb and lc, in which the superconductor effectively forms part of an external reservoir and 
is not simply part of the scattering region. Figure Id indicates a third class of N-SS'-N structures, involving two (or 
more) separate superconductors S and S', with respective order parameter phases <p, </>' . As noted above, transport 
properties of such Andreev interferometers are periodic functions of the phase difference <f> — <f>' . Figure le shows a 
fourth class of S-N-S' structures, in which two superconducting reservoirs are connected to a normal scattering region. 
In this case, the structure forms a Josephson junction and in contrast with all other structures shown in figure 1, in the 
linear response limit, the dc conductance measured between the superconducting reservoirs is identically zero. In this 
case, the relevant dc quantity is the current-phase relation and the associated critical current. Clearly one could also 
measure the Josephson current between the two superconductors in the N-SS'-N structure of figure lc and therefore 
the above classification is intended to label the measurement being made, rather than the device being measured. 
Structures of the form le will not be discussed. 

For the most part, the theoretical descriptions discussed below have finessed problems of self-consistency, by com- 
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puting measured quantities as a function of the superconducting order parameter pairing field f aa '(r) — (Vv(i:)Vv(z:)) 
induced by making contact with a piece of superconductor. As an example, figure la shows a normal mesoscopic 
scattering region in contact with a superconducting island which plays the role of an externally controllable source of 
/ctct'(e)) m much the same way that the coils of a magnet are an external source of magnetic field. The coils are not 
of primary interest and in many cases, neither is the superconductor. It is assumed that parameters characterizing 
the superconductor are given and the key question is how does superconductivity influence transport through the 
scattering region. Of course once the influence of superconductivity is understood, transport properties can be used 
to probe the symmetry and spatial stucture of the order parameter, as suggested by Cook et al. [1996]. Furthermore, 
in the presence of large currents which modify the order parameter, a complete self-consistent treatment is necessary, 
as described for example in [Bruder 1990, Hara et al 1993, Barash et al 1995, Canizares and Sols 1995, Chang et al 
1995, Gyorffy 1995, Martin and Lambert 1995 and 1996, Riedel et al 1996]. 



II. THE MULTIPLE SCATTERING APPROACH TO DC TRANSPORT IN SUPERCONDUCTING 

HYBRIDS. 

A. Fundamental current-voltage relations. 

In this section we review the multi-channel current-voltage relations for a disordered phase-coherent scatterer 
connected to normal reservoirs, obtained for two normal probes by [Lambert 1991] and extended to multi-probes by 
[Lambert Hui and Robinson, 1993]. The coefficients (denoted Aij and below) are directly accessible experimentally 
and can be combined to yield a variety of transport coefficients, the simplest of which is the electrical conductance 
of a normal-superconducting interface. To avoid time-dependent order parameter phases varying at the Josephson 
frequency, which would render a time-independent scattering approach invalid, these are derived under the condition 
that all superconductors share a common condensate chemical potential fi. The derivation of the fundamental current- 
voltage relation presented in [Lambert 1991] follows closely the multi-channel scattering theory developed during the 
1980s for non-superconducting mesoscopic structures [Biittiker 1986, Buot 1993]. In the normal state, this approach 
yields, for example, the multi-channel Landaucr formula [Landauer 1970] for the electrical conductance 

G = (2e 2 /h)T , (2) 

where To is the transmission coefficient of the structure. Historically the above formula was not accepted without 
a great deal of debate and as we shall see in the following section, contradicts the corresponding expression used by 
practitioners of quasi-classical theories, where the alternative expression 

G = (2e 2 /h)(T /R ) (3) 

is employed, with (in one-dimension) Rq the reflection coefficient. In fact the above two expressions refer to a 
two-probe and to a four-probe measurements, respectively, and the crucial lesson from the debate surrounding these 
equations is that transport coefficients such as the electrical conductance are secondary quantities. More fundamental 
are the current-voltage relations describing a given mesoscopic structure. 

Equations (|J) and ([}]) arc not valid in the presence of Andreev scattering, because charge transport and quasi- 
particlc diffusion are no longer equivalent. For example when a quasi-particle Andreev reflects at an N-S interface, 
the energy and probability density of the excitation are reflected back into the normal conductor, whereas a charge 
of 2e is injected into the superconductor. Thus charge flows into the superconductor, even though the excitation 
does not and as a consequence, a current-voltage relation should be used, which takes into account this charge-energy 
separation. For a scattering region connected to L normal reservoirs, labelled i = 1, 2, ... L, it is convenient to write 
this in the form 

N 

* = (4) 

3=1 

where Ii is the current flowing from reservoir i and the coefficients Ay will be discussed in detail below. In the 
linear-response limit, this reduces to 

L 

It = ^a.jivj ■- v), (5) 
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The above expressions describe reservoirs at voltages Vi,i — 1, 2, . . . , L, connected to a scattering region containing 
one or more superconductors with a common condensate chemical fi and relates the current Ii from reservoir i 
to the voltage differences (vj — v), where v — fi/e. The L = 2 formula describes a wide variety of experimental 
measurements and underpins many subsequent theoretical descriptions of disordered N-S interfaces and inhomogeneous 
structures. For this reason, after discussing the relationship between the coefficients Aij , and the scattering matrix, 
we shall examine the two-probe formula in some detail and illustrate its application to some generic experimental 
measurements. 



B. Relationship between the generalised conductance matrix and the s-matrix. 

The L — 2 analysis of [Lambert 1991] is based on the observation that in the absence of inelastic scattering, dc 
transport is determined by the quantum mechanical scattering matrix s(E,H), which yields scattering properties at 
energy E, of a phase-coherent structure described by a Hamiltonian H. If the structure is connected to external 
reservoirs by open scattering channels labelled by quantum numbers n, then this has matrix elements of the form 
Sn,n'{E, H). The squared modulus of Sn,n' {E, H ) is the outgoing flux of quasi-particles along channel n, arising from 
a unit incident flux along channel n'. Adopting the notation of [Lambert, Hui and Robinson 1993], we consider 
channels belonging to current-carrying leads, with quasi-particles labelled by a discrete quantum number a (a = +1 
for particles, —1 for holes) and therefore write n — (I, a), where / labels all other quantum numbers associated with 
the leads. With this notation, the scattering matrix elements s n . n '{E, H) = s^f{E, H) satisfy the unitarity condition 
s'(E,H) = s~ 1 (E, H), the time-reversibility condition s t (E,H) = s(E,H*) and if E is measured relative to the 
condensate chemical potential ^ = ev, the particle-hole symmetry relation s?',,(E,H) — a/3[s7u'~ (— E, H)]* . (For 
convenience we adopt the convention of including appropriate ratios of channel group velocities in the definition of s 
to yield a unitary scattering matrix.) 

For a scatterer connected to external reservoirs by L crystalline, normal leads, labelled i = 1, 2, . . . , L, it is convenient 
to write I = (i, a), V — (j, 6), where a(b) is a channel belonging to lead With this notation, the quantities entering 
the current-voltage relation are of the form 

P?/{E,H) = £ \s^ amb) (E,H)\ 2 = Tr [ S ^(E,H){sf(E,H)}^}, (6) 

a,b 

which is an expression for the coefficient for reflection (i = j) or transmission (i ^ j) of a quasi-particle of type (3 in 
lead j to a quasi-particle of type a in lead i. For a ^ (3, P?f(E,H) is an Andreev scattering coefficient, while for 
a = P, it is a normal scattering coefficient. Since unitarity yields 

0bj aaj 

where i and j sum over all leads containing open channels of energy E, this satisfies 

j2Pt/(E,H)=N?(E), and £ pf{E, H) = N?(E), (8) 

0j ai 

where N?(E) is the number of open channels for a- type quasi-particles of energy E in lead i, satisfying N^~(E) — 
N~(—E). Similarly particle- hole symmetry yields 

P?/(E,H)=P^>-P(-E,H) (9) 

and time reversal symmetry 

P^{E,H) = P^ a {E,H*). (10) 

Having introduced the scattering coefficients P^ (E, H), the coefficients Aij of the fundamental formula (^) are 
given by 

/>oo 

^ = (2e/fc)5>) / dE{8 ij N^E)f?{E)-Y J P^{E,H)ff{E)}, (U) 

„ JO a 
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with ff{E) = {exp[(E — a(evj — fi))/kbT] + the distribution of incoming a-type quasi-particles from lead j. [] 
Equation (||) yields the current-voltage characteristics of a given structure at finite voltages, provided all scattering 
coefficients are computed in the presence of a self-consistently determined order parameter and self-consistent values 
of all other scattering potentials. At finite temperature, but zero voltage, where vi — v — > 0, it reduces to equation 
(§, with a,ij given by 

a ij = (2e» J°° dE[-^l][N+(E)6 ij -P±+(E,H)+P^(E,H)], (12) 

where f(E) is the Fermi function and equation (|^) has been used. At finite voltages, but zero-temperature, it reduces 
to 

A V} = (2e/h) / dE [S ZJ N+(E) + Pr+(E, H) - P++(E, H)}. (13) 

Jo 

Finally at finite voltages, the differential of equation (||) with respect to Vj (with fi and all other potentials held 
constant) yields 

dli/dvj = aij, (14) 

where at finite temperature, 



aij = (2e 2 /ME( a )/°° dElSiiNrm-a^^p-] -^Pf^H)^-^!]}. (15) 
and at zero temperature, 



oy - [2e 2 //i][^A+(^) + Pr+{Ej,H) - P± + {Ej,H)]. (16) 

where where Ei = evi — fi. 

It is worth noting that replacing E by —E and utilizing the particle-hole symmetry relation (^|) allows equation 
( |l2l ) to be rewritten in the form 

aij = (2e 2 /h) |" dE[-^l}[Nr(E)S ij -P i :-(E,H) + P+-(E,H)], (17) 

which demonstrates that particles and hole are treated on an equal footing in equations ( |l2| ) and (p"7|). Furthermore 
in view of the symmetries (^), ( pp| ) the reciprocity relation aij(H) = aji(H*) is satisfied. 

C. Two probe formulae in more detail. 

While the above notation is convenient for arbitrary L, it perhaps obscures the simplicity of the final result and 
therefore in the literature, several alternative notations have been employed. For the case of L = 2 normal probes, 
where the scattering matrix has the structure 

(is) 

it is convenient to write s^(E,H) = r a/3 (E), s^{E,H) = r' af3 (E), s%?(E,H) = t afj (E) and s^{E,H) = t' af3 (E). 
With this notation, the sub-matrices r,t,r',t' have the form 

•-<*> = few 11$) • < 19) 



1 It is perhaps worth noting that in [Lambert 1991], the following notation is employed N P (E) = N± , R pp = 
P++(E,H)/N P (E), R hp = P{i+(E,H)/N„(E) etc. 
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The matrix r Q ^ (r' a p) is a matrix of amplitudes describing the reflection of /3-type quasi-particles from reservoir 1 (2) 
into a-type quasi-particles travelling back into reservoir 1 (2). Similarly t a p (t' aj3 ) is a matrix of amplitudes describing 
the transmission of /3-type quasi-particles from reservoir 1 (2) into a- type quasi-particles of reservoir 2 (1). 
Using these sub-matrices, the L = 2 current-voltage relation can be written 



where 



with 



and 
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/Trace {r ++ (E)rl + (E)}\ 
Trace {t ++ (E)tl + (E)} 
Trace {r_+ 

V Trace {t-+(E)tl + (E)} J 



/Trace {r' ++ (E)r'l + (E)}\ 
Trace {t' ++ (E)t'l + {E)} 
Trace {r'_ + (E)r'l + (E)} 

V Trace {t'_ + (E)t'l + (E)} J 



(22) 



(23) 



Similarly the zero temperature differential conductance (|16| ) becomes 

Oij = [2e 2 /h][5^N+{E t ) + R a {E 3 ) - R^E,)]. 



(24) 



D. Applications of the two-probe conductance matrix. 

As noted in [Lambert 1991], the two-probe current- voltage relation (|2^) can be used to derive generalisations of both 
and (|^). To further illustrate the versatility of the two- probe theory, we now apply it to some typical experimental 
measurements. For convenience in what follows, we set 2e 2 /h equal to unity. 

Example 1. The normal limit. 

Consider the structure of figure la, where a scattering region containing superconducting inclusions is connected 
to normal reservoirs at potentials v\, v%- In the normal limit, the condensate potential v must disappear from the 
fundamental current- voltage relation (^) . As shown below, this occurs, because in the absence of Andreev scattering, 
unitarity of the scattering matrix implies that an = — a\2 = —021 = 022 = Tq. Hence equation ( po|) reduces to 

h = -h = T Q {v x - v 2 ), (25) 

which is simply the Landauer formula (||). 

Example 2. Experiments where /12 = A*- 

Figure lb shows an experiment in which the superconductor and reservoir 2 are held at the same potential. In this 
case, equation (|2^) yields 

t-^t = a u = f° dE - ^1[N+(E) - R (E) + R a (E)}. (26) 
(ui - v) J_ oc dE 

This experimental configuration is of the type used in tunneling experiments, aimed at probing the proximity effect 
in the vicinity of an N-S boundary [Gueron et al 1996]. This result describes any of the structures shown in figure 
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1, provided fi 2 = M and is a generalization to disordered and inhomogeneous structures of the boundary conductance 
formula derived by [ Blonder, Tinkham, Klapwijk 1982]. 

Example 3. Experiments where ai2 = a.21 = 0. 

An example of such an experiment is shown in figure lc, where under sub-gap conditions, the presence of a long 
superconductor (of length greater than the superconducting coherence length) prevents the transmission of quasi- 
particles from reservoir 1 to reservoir 2 and vice versa. In this case, combining the unitarity condition N^(E) = 



Ro(E) + R a (E) with equation (|20|) yields 



h = an = r dE (-^W- ) 12R*(E)}- 



In common with the current-voltage relation from which it derives, equation ( p7| ) is valid in the presence of disorder 
and inhomogeneities and in the presence of an arbitrary number of superconducting inclusions of arbitrary geometry. 

Example 4. Experiments where Ii = I 2 = I. 

Such a situation is illustrated in figure la. In this case, inverting equation (|2^) yields 

Vi-v\ _ 1 / a 22 -a 12 \ ( I \ 

v 2 -v)-d[-a 21 ail )[-l)> (28) 
where d = a\\a 22 — 012021- Hence the two-probe conductance G = I /{v\ — V2) takes the form 

G= . (29) 

on + a 2 2 + 012 + 0,21 

As an example of this formula, we note that in the zero-temperature limit, where all quantities are evaluated at 
zero energy, equation (|2|) can be written [Lambert 1993, Lambert, Hui and Robinson 1993] 

For a symmetric scatterer, where primed quantities equal unprimcd quantities, this reduces to G = Tq + R ai whereas 
in the absence of transmission between the reservoirs, the resistance G~ x reduces to a sum of two resistances G _1 = 
(l/2R a ) + (l/2R' a ). It should be noted that a combination of particle-hole symmetry and unitarity yield at E = 0, 
To + T a = Tq + T^ and therefore equation ( p0| ) is symmetric under an interchange of primed and unprimed coefficients. 

Example 5. Experiments where I2 = 0. 

As a final example, consider the experiment sketched in figure Id where reservoir 2 acts as a voltage probe, with 
I2 = 0. In this case equation (^0|) yields for the ratio of the voltages 

02 -V) _ 021 

7 r — > 

(Vi- V) a 2 2 

where from equation (|2l]), the coefficient a 22 is positive. In contrast, the coefficient a 2 i is necessarily negative for 
a normal system, but in the presence of Andreev scattering can have arbitrary sign. Hence superconductivity can 
induce voltage sign-reversals which are not present in the normal limit. This feature was first predicted within the 
context of negative four-probe conductances [Allsopp et al 1994] and has been confirmed in recent experiments by the 
Groningen group [Hartog et al 1996]. 



E. The Bogoliubov - de Gennes equation. 

The above formulae relate measurable quantities to scattering matrix elements and therefore to end this section 
we briefly introduce the Bogoliubov - de Gennes equation [de Gennes 1989], which forms a basis for computing the 
scattering matrix s. The Bogoliubov - de Gennes equation arises during the diagonalization of the mean-field BCS 
Hamiltonian, which for a non-magnetic, spin-singlet superconductor takes the form 
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where Eq is a constant, Vv(n) and t/4-fc) are field operators, destroying and creating electrons of spin a at position r 
and 

H(rrl) _(5(r-r>)H (r) A(r,r') \ 

^II,Z1J-^ A*(r,r') S(r-r')H *(r)J ■ [66) 

(u (r)\ 

To each positive eigenvalue E n of if, with eigenvector *„(r) = W I satisfying 



d^H(r^)[:r^)=E n (: r ^), (34) 



/ "n 

v n (r') J ~™ \ v n (r 

f-v*(r)\ 

there exists a corresponding negative eigenvalue — E n with eigenvector vE'_ n (r) = *} \ )■ Consequently H e ff is 

\ u n\L) J 

diagonalized by the transformation 

V'tfc) \ _ V- f MD ~v*(r) \ ( 7«T \ (o^ 

4(t))~^M(r) <(L))W n J> (35) 

where to avoid overcounting, only one of }&_ n (r) or }&__ n (r) is included in the sum over n. 

It should be noted that whereas the Bogoliubov - de Gennes equation (Q) refers to a model in which the variable 
r varies continuously, in a tight-binding model, often used in numerical simulations, the corresponding Bogoliubov - 
de Gennes equation is 

Efr= -e^+ 7 *^-0 l+5 + A*^ [ ' 

where ipi indicates the particle (hole) wavefunction on site i and i + S labels a neighbour of i. 

As discussed in [Hui and Lambert 1990] and [Lambert, Hui and Robinson 1993], for a scattering region connected 
to two crystalline normal leads, the Bogoliubov-de Gennes equation may be solved by means of a transfer matrix 
method, which yields a transfer matrix T satisfying 

t(*) , (37) 



r j \o / 

where O (I) refer to vectors of outgoing (incoming) plane- wave amplitudes on the left and O' , (/') to corresponding 
amplitudes on the right, each plane-wave being divided by the square root of its longitudinal group velocity to ensure 
unitarity of s. Whereas T connects plane wave amplitudes in the left lead to amplitudes in the right lead, the s-matrix 
connects incoming amplitudes to outgoing amplitudes and satisfies 



0\_ I 
O' j I I' 



(38) 



Once T is known, the s-matrix can be constructed. Indeed if s is written as 

(39) 



r t' 
t r' 



then T has the form 



I'll T„ \ _ I ((>)-■ r' (*')-' 



t =<k; aJ-i-OT-r (?>-■ 1 • (40) 



from which the following inverse relation is obtained, 



Tp — lrri rp — 1 

22 J 21 J 2 2 



(41) 



An alternative method of evaluating coefficients in the above current-voltage relations is provided by the recursive 
Green's function method, which uses Gaussian elimination to compute the Green's function on sites located at 
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the surface of the interface between external normal leads and the scatterer. Given the surface Green's function, 
scattering coefficients can be obtained from generalised Fisher-Lee relations [Fisher and Lee 1982], derived by [Takane 
and Ebisawa 1992(a)] and later rederived in [Lambert 1993, Lambert, Hui and Robinson 1993]. This recursive 
technique is identical to the "decimation" method employed by [Lambert and Hui 1990] and is essentially an efficient 
implementation of Gaussian elimination. 

A third method of evaluating the coefficient R a was derived by [Beenakker 1992] for the case where there is perfect 
Andreev reflection at the boundary of a clean superconductor. The resulting formula expresses R a in terms of 
scattering properties of the normal state, as discussed in section V below and facilitates the application of random 
matrix theory to N-S structures [Beenakker 1997]. 



III. ZAITSEV'S BOUNDARY CONDITIONS AND MULTIPLE SCATTERING TECHNIQUES. 

In the previous section, the multiple scattering approach to dc transport was reviewed and in the following section 
the method of quasi-classical Green's function will be discussed. Whereas the multiple scattering method is a stand- 
alone approach, the quasi-classical technique must be supplemented by boundary conditions describing quasi-particlc 
scattering at the interface between different metals. These boundary conditions are provided by scattering theory 
and therefore provide a bridge between the two approaches. In the literature, quasi-classical boundary conditions 
have not been discussed using the language of modern scattering theory and as a consequence are restricted to planar 
surfaces with translational invariance in the transverse direction. In this section we offer a rederivation of Zaitsev's 
boundary conditions [Zaitsev 1984], which not only fills this gap, but also yields a more general condition applicable 
to non-planar interfaces. 

As noted in section 2, the starting point for a multiple scattering description is the scattering matrix s, with matrix 
elements s nn < connecting incoming to outgoing scattering channels, and satisfying equation ([38]). 

Clearly s(E,H) is a functional of all physical potentials entering the Hamiltonian H, as well as a function of E. 
Since H is Hermitian, quasi-particle probability is conserved, which yields 

s- 1 (E,H) = s^E } H) (42) 



' !1 ' T(E, 11)1(1 _°j r(£,fl). (43) 



and 

(o —j ; " \ o -i j 

Furthermore time reversal symmetry yields 

s*(E,H*) = s-\E,H) (44) 

and 



T(E,H)=(° 1 ^ T*(E,H*) 



(45) 



Hence 



and 



s(E,H*) = s l {E,H) (46) 



T- 1 \E,H*)=( l ljT(E,H)(J J). (47) 

It should be noted that the the unitary matrix s in equation (pq ) and the matrix T in equation (j3?]) connects open 
channels to open channels only. If evanescent states which decay at large distances from the scatterer are included 
in the outgoing states on the left-hand-side of (|3^) and states which grow at large distances are included on the 
right, then the corresponding scattering matrix s and transfer matrix T will not satisfy equations (E2) and (f43|). The 
matrices s and T should then be constructed by eliminating closed channels from s and T. 

Using the above definitions, the boundary conditions derived by Zaitsev are readily extended to the case of multiple 
scattering channels. The problem to be solved is that of connecting the Green's function Gj in a region j to the Green's 
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function(s) in all other regions connected to the scatterer. Although the following analysis is easily generalized to 
many leads attached to a scattering region, for simplicity of notation we restrict the discussion to only two leads. 

For the purpose of this section, it is convenient to write the channel index n in the form n — j, <r,p, where j = 1,2 
labels a lead, p labels all tranverse quantum numbers, a = +1 for a right-going (or right-decaying) wave and a = — 1 
for a left-going (or left-decaying wave). If denotes a position in a crystalline lead j of constant cross-section, then 
a plane wave of unit flux incident along channel n can be written 

ULj) = Xn(fi) K(S))- 1/2 exp( l fcf CT (i?)z J ), (48) 

where Xn{p ■ ) is the transverse mode associated with channel n and (v n (E))~ 1 / 2 exp(ik?(E)zj) is a plane wave of unit 
flux, with longitudinal wavevector k? a (E). 

Similarly for ^ r,-, and r,*, r_j located within the leads, the Green's function has the form 



V«(i,<r, P )(i>o>>o(£) 



i expi[kl (T (E)z i + kl al (E)z j 

where £ is equal to the sign of Q&i — Zj). 

For r { j^z r_. and a fixed value of S, the Green's function is simply a wavefunction and therefore the matrix A % ^ a , with 

matrix elements A 1 ^ a '(p,p' satisfy relations of the form of equation ([37]). First consider the form of the Green's 
function when Zj < Zi and j — 1. When viewed as a function of Zj, for Zj < 0, one has 

/All All \ (All All \ 

T(E,H)( A A A 7_f_)={X %Z) (50) 

On the other hand, when viewed as a function of Zj, using the conjugate equation for the Green's function (which 
involves the time reversed Hamiltonian, one obtains for i = 2 and Zi > 0, 

422 422 \ / 421 421 \ 

A? a™ + ) =T{E > H * ) [a¥ A^ + ) (51) 

Hence after eliminating the off-diagonal terms, A 2 },, etc., we obtain the general boundary condition relating the 
Green's functions on the left of the scatterer, to the Green's functions on the right: 

/A 11 A 11 \ ( A 22 A 22 \ - 
T(E,H)(^l A^) = {a 2 \ ^Tjm-ff*)- 1 - (52) 

Equation ( [52] ) is a generalisation of Zaitsev's boundary condition to the case of a non-planar barrier, which in 
general may contain impurities and break time-reversal symmetry. ^| 

Since the sum on the right-hand-side of equation ([l9|) includes evanescent channels, the general boundary condition 
involves the matrix T. However at large distances from the scatterer, evanescent channels can be ignored and therefore 
T can be replaced by T. In view of equation (f47|), if time-reversal symmetry is present, the boundary condition 
simplifies to 

/-A 11 A 11 \ ( -A 22 A 22 \ 
T(E,H)^tT A%_) = {-%1 aK) T{E ' H) (53) 

In the work of Zaitsev, where p represents the transverse momentum and all matrices are diagonal, equation (|5 
is satisfied for each separate p and in the notation of Zaitsev has the form 



A similar argument with Zj > Zi yields an identical result and therefore the boundary condition is independent of the choice 
of E. 
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As an example, taking the trace of equation (|54j) yields |^] 

Si- Si =9% -9-, (55) 

which demonstrates that the anti-symmetric part of the quasi-classical Green's function is continuous across a scat- 
tering region More generally, equation ( |53| ) yields the corresponding relation 

Tr( A" + - A$_ ) = Tr(A 22 + - Af_ ) . (56) 

Similarly generalisations of all other boundary conditions can be expressed as traces over the matrices Af , . Indeed 
multiplying both sides of ( |53"|) by each of the matrices 



\ (T^ OW 0\ fO T^i 



T^ 1 ' I OJ' lTf 2 x or U 



(57) 



and taking the trace of the resulting four equations yields 

Tr[A n + - A 2 \] = -■R[2£ 1 T 2 iA£ f + T l2 T^ x A 2 3_] = Tr[rA^ + - r' A 22 _] (58) 

Tr[A£_ - Af_] = -Hrpr^Tia^ + T 21 T{ 1 1 Af + ] = Tr[rU n _ - r' ] A 22 + ] (59) 

TrL4 n + + A 2 ^] = TrhTf/Tn^V + T^Tf/AfJ = Tr^A^ + r 1 ' 1 A 22 + ] (60) 

Tr[A^_ + A 22 + ] = Tk[-T a l 1 r 2 ai4i 1 _ + TuT^ A 22 _] = Tr^Mi^ + r'^ 1 A 22 _] (61) 
Subtracting (|5q) from (E9I) yields, in view of (Ef 



Tr[rAW - r^A u _} = Tr[r'A 2 _ - A 22 + ] (62) 



and subtracting ( |60| ) from (61) yields 



Tr^A^ - r-M n _] = Tr^A 2 ^ - r'^AfJ (63) 
Adding (H) to (H) and @ to © yields 

Tr[A^_ + A n +] - TrLAf^ + A 22 ^ = Tr[r^+ + ^ A n _] - Tr[r'A 22 _ + r A A 22 + ] (64) 

and 

Tr[A%_ + A l \] + Ti[Af_ + A 22 + ] = Tr^A^ + r" 1 ^] + Tr^A 2 ^ + r* -1 A??_]. (65) 

For the translationally invariant case considered by Zaitsev, writing t = t' = \t\expi9 r = |r|expi</>, r' = 
— \r\ expi(26> — 0), A" CT = Ql and A"_ CT — g£, one obtains from either of equations (^3) and (^3) 

<?a = & (66) 



• ! # = M(#+#), (67) 



from m§ and (B5J) 



= ^(Gl-Gfh (68) 



3 More precisely, one has to multiply eq. (|54|) by T 1 and exploit the ciclyc property of the trace. 
4 As it will become clear in the next section, this corresponds to the conservation of the current 
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and from ( |55| ) 



al = al (69) 

where 

Ola = [Q\ exp i<j> ± 0i exp -^]/2, (70) 



= [0+ exp - 20) ± G 2 _ exp i(20 - <f>)}/2 (71) 

and 

£ ia =$_±£_]/2. (72) 

Equations ( |66| ) and ([39]), which are a limiting case of the more general relation ( (53| ) are identical to the boundary 
conditions of Zaitsev and will be discussed in the next section. 



IV. QUASI-CLASSICAL GREEN'S FUNCTION APPROACH 

A. Quasi-classical equations 

The quasi-classical approach to the theory of superconductivity, initiated by Eilenberger [1968] and Larkin and 
Ovchinnikov [1968] and developed by several authors [Usadel 1970, Eliashberg 1971, Larkin and Ovchinnikov 
1973,1975] has been largely used to analyse transport phenomena in dirty hybrid systems. This approach can be 
used to study thermodynamic and kinetic properties of superconductors, whose dimensions significantly exceed the 
Fermi wavelength Xp = 27r/fcp. For the purpose of reviewing recent theoretical work, we briefly summarize the 
quasi-classical Green's function approach, following mainly [Larkin et Ovchinnikov 1986]. There exist in the literature 
several excellent reviews on the subject including [Rammer and Smith 1986] and we refer the reader to these for a 
more detailed exposition. The derivation of the equation for quasi-classical Green's function starts from the Dyson 
equation for the matrix Green's function G in the non equilibrium Keldysh formalism [Keldysh 1964)]. 

(Go" 1 - t)G = 1 (73) 

where 



Following standard notation, we indicate with a "hat" matrices in Nambu space and define retarded, advanced and 
Keldysh Green's functions 

G*(l, 2) = 9(h - ta)(G-+(l, 2) - G+" (1, 2)) (75) 
G A (1,2) = -6(t 2 -t 1 )(G- + (l,2) - G+-(l,2)) (76) 
6(1,2) =G-+(1,2) + G+-(1,2). (77) 



where 



■r-+ = ( <Vr(i)V»J(2)> <^Wj( 2 )> i (m 

" 1 -<^(1)V>{(2)> -< ^(1)^(2) > 1 



;r+- = _( <^{(2)Vt(i)> <V;(2M(i)> \ r7ql 

l - < v4 (2)^(1) > -< ^(2)vl(i) > r 1 ] 
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In these expressions, 1 = (ri,ii) and 2 = (^,£2), which allow equation ( |73| ) to be written in the less compact form 

d2(&o 1 (l, 2) - S(l, 2))G(2, 3) = 5(1 - 3). 

The equation conjugate to equation ( |73| ) is 

G(Gq 1 - t) = 1 (80) 
and taking the difference between equations ((73j) and ( |so"| ) yields 

[Go 1 -£,G]=0. (81) 
This equation is simplified by going to the center-of-mass and relative coordinates (R, T) and (r, t) defined as 

r 1)2 = R±r/2, t 1 , 2 =T±t/2, 
Fourier transforming with respect to r and £, and introducing the quasi-classical Green's function defined by 

p CO 

g(H,T ] p,e) = - d£G(R,T;p,e) (82) 



where £ = p 2 /2m — fi is the energy measured with respect to the Fermi level and p is the unit vector in the direction 
of the momentum p. On the assumption that the self-energy depends weakly on the energy £, one sets £ = in E, 
which yields 

9t{t z ,§} + v F p ■ d-R,g - it [f z ,g] + i [E, g] = 0. (83) 

where the curly brackets denote the anticommutator. This is the equation for the quasiclassical Green's function 
first derived by Eilenberger [1968]. f z is a block-diagonal matrix with diagonal block entries a z . In deriving equation 



(30), the non-homogenous term present in the Dyson equation (73) was eliminated and as a consequence, equation 
(33) determines g only up to a constant. As discussed in [Shelankov 1985, Zaitsev 1984], a useful normalization 
condition for g is 

99=1- (84) 
In terms of the quasiclassical Green's function, the physical current is given by 

1 f 00 

j(R, T) = —N vf / de < [iTr{a z g(e,p, R, T)) >, (85) 

where iVo = m 2 vp/2i:' 2 is the free single-particle density of states per spin, a z the usual Pauli matrix and vf 
the Fermi velocity. Here < ... > indicates the average over the angle 8 formed by p and the direction of R, i.e., 
(1/2) dcos(#)(...). It should be noted that, because of the angular average, only the antisymmetric part of g 
enters the expression for the current. 

In the dirty limit, in the case of an isotropic scattering potential with elastic scattering time t, the effect of 
non-magnetic impurities can be described by the self-energy 

S = < 9 > , (86) 
It 

and equation ( |83| ) can be considerably simplified [Usadel 1970]. In this case one expands g in spherical harmonics 
keeping only the s- and p-wave terms 

ff(cos(0)) =g + cos(f?)3! (87) 

with go and g\ not depending on cos{9) and cos(8)g\ <ti go- By inserting equation ( [8^ ) in equation ([83]), g\ is expressed 
in terms of go 
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91 = -Igodngo 



(88) 



and for go one obtains a diffusion-like equation 

Dd R g d R g + ie [f z ,g Q ] - d T {f z ,g } = (89) 

with D = vfIt/3 the diffusion coefficient and I = vpr the mean free path. 

Eqs.([B3l) and ([89] ) are valid for bulk systems and in the presence of boundaries, must be supplemented by boundary 
conditions which connect the Green functions evaluated in different regions. These boundary conditions have turned 
out to be crucial to the recent development of transport theory in hybrid structures. For this reason, we present a 
brief discussion of these in the following subsection. 



B. Boundary conditions for quasiclassical Green functions 

In this subsection, we discuss the boundary conditions for the quasiclassical Green's function [Zaitsev 1984], given 
by equations (j6(^) to (^9|) of section III. The antisymmetric functions g a — g\~ g 2 and G a = Gl = G 2 are continuous 
across the boundary, while the symmetric ones g\ and Q\ experience a jump determined by the transparency of the 
barrier. The size of the jump vanishes for perfectly transmitting interfaces. 

The function g is the quasi-classical Green's function satisfying equation whereas the Green's function Q 

arises from reflected waves at the boundary. It is therefore necessary to derive an equation of motion for Q to be used 
together with equation (jS3[). We will not show this derivation here, because in actual calculations, one only uses the 
Green's function g. Zaitsev has shown, that the final result is of the form 

g\h)g\h) = (-lysgnik^gHh), (90) 

which must be used in conjunction with the normalization condition 

.'/'.'/' = I. (91) 

The boundary conditions for the antisymmetric components are already decoupled with respect to the g a and Ga- 
To decouple the symmetric components, we express eqs.(|9l^-|9~l)) in terms of the symmetric and antisymmetric parts 
to yield 



9lGl + 9aGa = (-lYGa 
9 i s Ga + 9aG i s = {-l) i G i s 



(92) 



and 

glgl + 9a9a = i, g\g a + g a gl = 0. (93) 

By manipulating equation (|92] ) one obtains 

glG] + &G\ = Ug\G\ - g 2 s Gs) (94) 



which is the extra condition to be used together with equations (|67[]68|) . By means of equations Q67j-|SS|) one expresses 
G 2 and Gl in terms of g\ and g 2 and substitute in equation (]94|) so that the final boundary condition reads 

g a [R(l - g a g a ) + (T/4)(g s 1 - g 2 ) 2 ] = (T/4)(g 2 jl ~ glgl) (95) 

where R = |r| 2 , T = 1 — R are the reflection and transmission coefficient of the barrier. 

As a simple example, we consider the normal case, when the Green's function is a two-by-two matrix 



i g 
o -1 



(96) 



with g = 2/, / being the usual distribution function entering the Boltzman kinetic equation. By observing that 

[g ls , g 2s ] = 4 (° n fls ~ hs ) ; ( 9l - g 2 ) 2 =6; g 2 a = 6 (97) 
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the boundary condition (|95j) assumes the form 

f a = (T/2R)(f ls -f 2s ). (98) 
It is worth noticing that this result could have been obtained directly from simple counting arguments of the form 
h{k) = Tf 2 (k) + Rh(-k), f 2 (-k) = Tfxi-k) + Rf 2 (k). (99) 



In view of equation (|35j) for the current, equation ( J98| ) yields for the conductance of a tunnel junction in the normal 

case, 

G T = e 2 N v F < fxT/R > (100) 

where \i = 6. In the general case, equation (Bq) can be simplified in the low transparency limit. If T <C 1, then 
g a ~ T and the boundary condition equation (|95|) reduces to 

g a =(T/4R)[g ls ,g 2s ]. (101) 

In the dirty limit, g a = —^lg s djig s and equations ( |S8f]8"9| ) yield the boundary condition [Kuprianov and Lukichev 
1988] 

lgi s d R g ls = lg2sd R g2 S , (102) 



which after multiplying equation (101) by /i and taking the angular average yields 



3 

Igisdngis = - < pTjR > [g 2s ,<?i s ]- (103) 



By defining a conserved "super" current / (cf. equation (p5|) ) one finally arrives at the following boundary condition 

(7 G'J 7 

I = -hsdnhs = — [fo s ,gu }, (104) 
e Ze 



where a = 2e 2 NoVpl/3 is the Drudc electrical conductivity. Equation (104) is the desired boundary condition to be 



used together with the diffusion equation (89) in the presence of boundaries. Eq.(104) is strictly valid in the case of 
small barrier transparency. Such restriction has been recently somehow relaxed by Lambert, Raimondi, Sweeney and 
Volkov [1997]. In the next subsection we show how to compute the current- voltage relation of the two fundamental 
elements of an hybrid system, namely a tunnel junction and a diffusive region. 

C. Quasi-classical theory at work. 

(i) A Tunnel junction. 

Consider the current through a tunnel junction as given by the commutator on the right-hand side of equation 



(104). The physical current is obtained from the Keldysh component of equation (104) 



[h,9i]k = 92 9i + h9\ " 9\h ~ 9\92 
where we have dropped the subscript "s". The normalization condition, gg = 1, allows us to choose 

.91,2 = 5^2/1,2 - fl,29l t 2 

where the matrix / can be taken to be diagonal [Larkin and Ovchinnikov 1975] 

fl.2 = /l.2<5o + f\p,O z . 

As a result, multiplying by a z and taking the trace, yields (cf. equation (pa)) 
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S 1 pOO 

3 = jf deTr(a z {fj a + fj b + ffl c + ffl d )) (105) 



where 



h = -[9m-9£)-{$-9i)9t] 

ic = [92(91^ - Oz9i) - (g?°z - z z g?)g A ] 

h = - [9x(92Vz - o z g A ) - (g R a z - a z g A )g A ] 
Due to the normalization condition g R ^ g R ( A ) = X, we have 

9 R ^=Z R W-a^9? W 



i=l 



where g R ( A *> = (iF R ^ sin(4>) , iF R ^ cos(4>), G R ^) and </> is the phase of the superconducting order parameter. 
Hence 



fi />oo 

3 = -f / de(Ij + I PI ) 



(106) 



where 



I.j = isinfa - <h) [f 2 (F 2 R - F A ){F R + F A ) + f°(F R + F A )(F R - F A )} 

and 

I PI = [(Gf - G A )(G R ~ G A ) + cos(^ - 4>2){F R + F A )(F R + F A )] (/* - /*). 
In equation (UOq), Ij is the Josephson current, while Ipj is sometimes referred to as quasi-particle and interference 



current. To clarify equation (106), consider the case where the two regions on the left and right of the barrier are 
at equilibrium and cj>i — <j) 2 , so that the Josephson current vanishes. If the distribution functions ff 2 have their 
equilibrium form, 

fi,2 = \{tanh{{e + eVi, 3 )/2T) - tanh((e - eV ll2 )/2T)) 



then the current through the junction becomes 

Gt 

where 



/oo 
de(ft - fi)M 12 (107) 
-OO 



M ia = |((Gf - G A )(G R - G A ) + (F R + F A )(F R + F A j). 

At T = 0, this reduces to 

I = G T M 12 U =0 {V 1 -V 2 ). (108) 



18 



Equation ( 108 ) shows that conductance of the tunnel junction is renormalized by a term M12 which depends on 
the amount of superconducting pairing on the two sides of the junction. In the limit of normal systems, = 1, 

and one recovers the conductance of the normal state (cf. equation (100)). 



(ii) A diffusive region. 



We derive now the equivalent of equation (107) for a diffusive region in the absence of inelastic processes. In this 
case the equation for the "hat" components of the Green's function (p9j) reads 



R(A) 



R(A) 



d R (g H d R g + gd R g A ) = 



Using g — g R f — fg A and the equation for g R ( A > , the equation for the Keldysh component becomes 



9r 



d*J - g H (d R f)g A - (g H d R g H )(d R f) - (d R /)(s^ R g 



0. 



The last equation yields /, once g R and g A have been determined. Choosing / = f°ao + f z <r z and g R = G R d z 
iF R a y , yields, after multiplying by a z and taking the trace, 



Or [(1 - G R G A - F R F A )d R f z ] = 0. 



(109) 



If we consider a diffusive region of length L, the distribution function depends only on the longitudinal coordinate 
x and f z {x) has the form 



where 



m(L) 



m{x) J dx> l^G R {x')G A {x')-F R {x')F A {x'Y 



(110) 



Using the formula for the physical current (pq) one obtains 



roc 1 



(in) 



To use equation (111) one has to solve the equation for g R . Since g R g R = 1, one can write G R = cosh(u) and 
F R — sinh(u), which yields 



F)d\u + 2iesinh(u) = 0. 



(112) 



Equation (112) is a non-linear equation, for which approximate analytical [Zaitsev 1990, Volkov et al. 1993, Volkov 
1994, Zaitsev 1994] and numerical [Zhou et al, 1995, Yip 1995, Nazarov and Stoof 1996, Golubov et al 1997] solutions 
have been obtained. For the purposes of the present discussion, we note that at zero energy, the solution may be 
obtained easily and takes the form 



, . u(L) - u(0) 
u(x) = v ' — —x + m(0). 



Since 



1 - G R G A - F R F A = 2cosh'{Re(u)), 



the current (111) becomes 



I _ a tanh(Re(u(L))) - tanh(Re(u(0))) 
~ L Re(u(L)) - Re(u(0)) ' 



(113) 
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where we have used the fact that for T — > 0, ^ de(f z (L) — / z (0)) — > V. The values and u(0) are 

determined by attaching the diffusive region to electrodes which are assumed to be in equilibrium. For a normal 
system u = 0, while for a superconducting one Re(u) = and Im(u) = — 1. As a result the resistance of a diffusive 
region does not depend on the nature of the electrodes. 

Nazarov recognized [Nazarov 1994] that the simplicity of the results of equations (108 113| ) can be generalized to 



an arbitrary circuit with normal and superconducting electrodes. In fact, the evaluation of the conductance can be 
reduced to the use of a set of simple rules, which we now state. A circuit is made of dispersive elements (diffusive 
regions and tunnel junctions) connected to one another and to the electrodes by nodes. To each circuit element is 
associated a "scalar" current, j, which is just the physical current and a "vector" current I. To each node is associated 
a "spectral" unit vector s. The vector current flowing through a dispersive element depends on the spectral vectors 
at the two ends. It is useful to imagine the spectral vector as the spherical coordinate in the northern hemisphere. 
The rules are: 

I) Conductance in the presence of superconducting electrodes is the same as in normal circuits, with renormalized 
tunnel junction conductances. The renormalization factor is given by the scalar product of the spectral vectors at the 
two sides of the junction. 

II) Spectral vector of a normal electrode is at the north pole (z), while that of a superconducting electrode lies on 
the equator, with the longitude given by the phase of the order parameter. 

III) Vector current is perpendicular to both the spectral vectors at the ends of a given element element. If <j> is the 
angle between the spectral vectors at the two ends of a given elements element, the magnitude of the vector current 
is I = Gd4> and / = GTsin(<p) for diffusive and tunnel junction elements, respectively. 

IV) Vector current is conserved at each node (Generalized Kirchoff rule). 
These rules will be applied to a N-I-S interface in the following section. 

V. THE N-S INTERFACE AND ZERO-BIAS ANOMALIES. 

In this section, we show how the multiple-scattering approach and quasi-classical theory can be used to describe 
the simplest example of phase-coherent structure, namely a N-S interface. 

A. BTK theory and the Andreev approximation 

The simplest problem of a one-dimensional N-I-S system was examined by [Blonder, Tinkham, and Klapwijk 1982]. 
This solution is also valid in higher dimensions, provided there is translational invariance in the direction perpendicular 
to the electronic motion. To see this, imagine a two-dimensional system with a finite width in y direction and infinite 
in the x direction. The motion along y is quantized, and for fixed k y , the problem becomes one-dimensional and the 
Bogolubov-de Gennes equations reduces to 

E1>(x) = {-^81 - /x)V(z) + A(xMx) 

where ~p = fj,— h 2 ky/2m is the effective chemical potential in a given channel. For a pairing potential with a step-like 
spatial variation of the form 

A(x) = A e(x), (115) 

the Bogolubov-de Gennes equations can be solved by matching wave functions. For x < the solution (ip, </>) for an 
incident electron from the left is 



and for x > 



ip L (x) = e tkx + r e 
<j) L (x) = r a e lqx 



ij) R (x) = t ue lkx +t a ue lqx 
(f> R (x)= t ve^ x + t a ve-^ x . 



(116) 



(117) 
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The coherence factors u and v correspond to a particle-like excitation of energy E and momentum k, while u and 
v correspond to a hole- like excitation at the same energy and with momentum q, where E = (H 2 /2m) fc 2 — ~p = 

-{h 2 /2m)q 2 +Jl and E 2 = [(?j 2 /2m)fc 2 - Jl} 2 + A 2 = [(h 2 /2m)q 2 +Jl} 2 + A 2 . Also 



u 2 {v 2 ) = $(l±y /W=5? /E) 

v 2 (u 2 )= 1(1 ± VE 2 - A 2 /E) [ j 

with sign(u/v) = sign(A/(E—ex)), sign(u/v) — sign(A/(E—e^)), where er = (h 2 /2m)k —~p and €q = (fi 2 /2m)q 2 + 
~p for £7 > A, while for £7 < A we have 



u/v = A/(E - ^ VA2 - E 2 ) 

u/v= A/(E + WA 2 - E 2 ). 1 J 

To solve for the scattering coefficients ro, r a , to, and t a one needs matching conditions at the interface x — 0. These 
are the usual continuity condition for (tp, </>) and their derivatives and read 

1 + r = t Q u + t a u 
r a = t v + t a v 

1 - r = ^t u - \t a u + (1 + r Q )U /tk (120) 
r a = jt v - lt a v + r a U /iq 



where the term in Uo in the last two equations of (120) allows for a delta-like potential U(x) — UoS(x) in order 
to reproduce the BTK model of an N-I-S interface. The system (120) completely determines the four scattering 
coefficients. The detailed solution is of course straightforward, but the resulting expressions look rather cumbersome, 
so it is convenient to resort to Andreev's approximation of setting k = q = k = q, which amounts to neglecting terms 
of the order A/E F . Writing z = 2mU /2kh 2 yields 



to = 


V 


(1 - iz 


) 


uv(l - 


t-*V 


vuz 2 


ta = 




ivz 




uv(l - 


-**)- 


vuz 21 


r a = 




!'?' 




uv(l 




vuz 2 


ro = 


(vu — 


uv)(z 


+ iz) 


uv(l - 


Vz 2 )- 


vuz 2 



(121) 



For E < A, where T = T a = 0, this yields 



E 2 + (l + 2z 2 ) 2 (A 2 - E 2 ) 
and in view of unitarity, Rq = 1 — R a . For E > A one obtains 

A 2 



A 2 

R a = w I Ts I o■.2^2^A2 ( 122 ) 



Ra — 



(123) 



{E + (1 + 2z 2 ) 2 \/E 2 — A 2 )) 2 
fl = 4z 2 (l + z 2 )(E 2 - A 2 ) 

° (£?+(! + 2z 2 )y/_ E 2 — A 2 )) 2 
T = 2£:(l + z 2 )(£^+ y/g 2 - A 2 ) 
° (£?+(! + 2z 2 ) y/g 2 - A 2 )) 2 

T = 2£;z 2 (£:- Ve 2 -A 2 ) 

a (E+(l + 2z 2 )(\J E 2 - A 2 )) 2 

The above treatment of a planar interface has been extended by a number of authors, including [Riedcl and Bagwell 
1993]. In one dimension and at zero temperature, the boundary conductance of equation ( p6| ) reduces to the BTK 
result 

G= (2e 2 /h)(l-Ro + R a ), (124) 
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which for E < A becomes G = (2e /h)2R a , while for E A, R a — > and one recovers the Landauer expression (2), 
G = (2e 2 /h)(l — Rq) = (2e 2 /h)To. It is interesting to compare the sub-gap conductance 



with that in the normal state 



G NS = (2e 2 /h) {i + 2z2)2 (125) 



G N = (2e 2 /h)T N = (2e 2 /h)-^. (126) 

1 + z z 



This allows one to express Gns m terms of the dimensionless normal state conductance Tn- 

2T 2 

G NS = (2e 2 /h) " (127) 
{l — In) 

which is a consequence of the fact that the superconducting electrode in this approximation only enters as a boundary 
condition. 



For a disordered N- region in contact with a superconductor the result ( 127 ) can be extended to higher dimen sion s 



by introducing the eigenvalues T n of the transm issio n matrix T — t pp t pp of the N-r egion and summing equation ( 127 ) 



over all n [Beenakker 1992], to yield equation (132) below. For z ^ 6 , equation (124) predicts that the low-energy 
conductance is suppressed by a factor (1 + z 2 )~ 2 , as compared to the z = case. Already for values of z ~ 3 this 
implies an almost vanishing sub-gap conductance. At the gap energy, the expression for R a achieves a peak value 
of unity within Andreev approximation, with R a decaying to zero for energies above the gap. This picture clearly 
resembles that of the tunneling approach, where the conductance is controlled by the supeconducting density of states, 
but is at variance with the experiment of [Kastalskii et al. 1991]. It was soon recognized that the low-energy peak 
observed by Kastalskii et al was due to Andreev scattering [van Wees, de Vries, Magnee, and Klapwijk, 1992], but 
it was not possible to explain it using the theory discussed above, because the interplay between scattering due to 
disorder in the normal region and Andreev scattering at the N-S interface is crucial. The first insight into the ZBA was 
based on a detailed description using quasi-classical Green function methods [Zaitsev 1990, Volkov and Klapwijk 1992, 
Volkov, Klapwijk and Zaitsev 1993]. Hekking and Nazarov [ 1993,1994], obtained similar results using a tunneling 
approach and Beenakker and co-workers [Beenakker et al 1994; Marmorkos et al 1993] confirmed these results using 
multiple scattering methods. 

To illustrate how these different methods explain the ZBA, the following sub-section briefly summarizes quasi- 
classical predictions for the zero-energy conductance and in sub-section C, the multiple scattering approach is dis- 
cussed. 



B. The quasiclassical approach: circuit theory. 

To apply the "circuit rules" of quasi-classical theory [Nazarov 1994] to a N-I-S system, we consider two resistances 
in series; the first is associated with the normal disordered region (N), while the second is associated with the 
tunnel junction (I). We assume that the normal region progressively widens in a macroscopic electrode and that the 
superconductor plays the role of the other electrode. In this circuit there are three spectral vectors. The vector 
associated with the normal electrode is sjy = z and that associated with the superconducting electrode is sg = x, 
both or which are fixed (cf . rule II of circuit theory) . The third spectral vector is located at the node connecting the 
diffusive element and the tunnel junction and lies in the x — z plane s = cos(#)x + sin(0)z. According to the circuit 
rules, one computes the conductance of this simple circuit in the usual way, but the resistance of the tunnel junction 
is renormalized (cf. rule I of circuit theory) 

Gns = -rjF, — , 1 ur< — ■ taw (^28) 
l/G D + 1/(Gt sm(0)) 

where Gu and Gt are the conductances of the normal diffusive region and of the tunnel junction in the normal state, 
respectively. The angle 9 is determined from the conservation of the vector current (cf. rules III- IV of circuit theory) 
and satisfies the equation 

G D 9 = cos(6»)G T . (129) 
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This equation has two simple limits: i) Gd/Gt » M « Gt/Gd, Gns ~ G t /Gd] ii) Gd/Gt < 1, 9« 7r/2, 
Gjvs ~ Gd + Gt- This shows that, by varying the parameter Gt/Gd, the conductance switches from a quadratic to 
a linear dependence on the tunnel junction conductance. A two-particle process behaves effectively as a one-particle 
process (Gns ~ Gt). The smallness of the probability for a two-particle transmission is compensated by the relative 
increase of disorder- induced scattering on the normal side. A detailed comparion between eq.( |l28| ) and a numerical 
S-matrix investigation can be found in Claughton, Raimondi and Lambert [1996]. 



C. The scattering approach: random matrix theory. 

To illustrate the multiple scattering approach to the ZBA, we now show how equation ( |128| ) may be derived using 
random matrix theory [Beenakker et al 1994; for a more detailed introduction to random matrix theory applied to 
quantum transport see also Stone, Mello, Muttalib and Pichard [1991] and the recent review by Beenakker 1997]. The 
starting point is the following scaling equation for the eigenvalue density p(X,s) —< 5Z„ =1 <5(A — X n ) > [P.A. Mello 
and J.-L. Pichard 1989], 

8 2 8 8 f°° 

gp(X,s) = -- — A(l + A)p(A, S )— j dX'p{X',s)\n\X-X'\ (130) 

where A„ is related to the transmission eigenvalue T n of the n-th transverse mode by 

A„ = (1 - T n )/T n (131) 

and the average < ... > is over the appropriate ensemble. Here s — L/l, where L is the size of the system and / the 
mean free path. N is the number of open channels in the normal region. The above equation yields the ensemble 
averaged eigenvalue density for a diffusive system and it is physically equivalent to the quasiclassical theory outlined 
in the previous section. Once p(X, s) is known, the conductance of the N-I-S system is given by (cf. eq.(127)) 

«™ = X g(23t? = X 2 X <"" (A ""(TW (132) 

With respect to the scaling variable s, the presence of a barrier at the N-S interface, acts as an "initial condition", 
instead of a "boundary condition" . In fact, at the interface, where s = 0, the eigenvalue density in a mode-independent 
approximation, is given by 

p(X, 0) = NS(X — (1 — T)/T), (133) 



where T is the transmission coefficient of the barrier. To solve the evolution equation (130) with the initial condition 
(133) one introduces the auxiliary function F(z,s) 

F(z,s)= HdX'?^, (134) 

Jo z — X 

which is analytic in the z-complex plane with a cut along the positive real axis. In terms of F, the evolution equation 
(130) becomes 

N—F+—z{l + z)F 2 = 0. (135) 
OS oz 



Writing z = sinh 2 (£) and 



yields 



U(Z,s) = S ^p±F(z(Z),s) (136) 



d d 

— U + U—U = 0, (137) 
os c>£ 
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which is Eulers equation of an ideal two-dimensional fluid. The real and imaginary parts of U are the cartesian 
components of the velocity field U — U x + iU y . The solution can then be obtained directly in terms of the initial 
condition Uq(£) = U(£,s), as 



U(Z,s) = U (Z-aU(t,a)). 



The initial condition (133) for the eigenvalue density translates for U in 

sinh(2£) 



2(cosh (£) - T _1 ) 



which in the limit of small transparency barrier, i.e., T <C 1 is approximated by (cf. eq.(lOl)) 

tfo(0=-§sinh(20- 



To evaluate the conductance, we note that in terms of the function U, the equation (132) becomes 



2e 



dU 



Gns = -r-N lim — 



(138) 



(139) 



(140) 



(141) 



By using (138) together with (140) and writing the real and immaginary part of U at £ = — i7r/4, we obtain the pair 
of equations 



U x 



T 

•— smh(2sU x ) sin(2sU y ) 



T 



U y — — cosh(2sZ7 x ) cos(2sU y ). 

The first equation requires U x = 0, and by defining 6 = 2sU y , we rewrite the second equation 

9 = sTcos{9). 



(142) 



(143) 



By recalling that — (2e 2 /h)(l/L) and Gt = (2e 2 /h)T are the conductances of the disordered normal region 
and of the tunnel barrier, we recover the equation for the conservation of the vector current ( |129| ) obtained within 
circuit theory. The expression for the conductance (141) may now be obtained by noticing that 



dU 



T cosh(2£ - 2sJ7) 



(144) 



and setting £ = — z7r/4. Insertion of equation(144) into equation (141), leads to an expression for the conductance 
identical to that given by circuit theory (cf. eq.(128)). 

The above ideas have been further developed in a number of treatments of Andreev scattering in chaotic and resonant 
stuctures, including [Altland and Zirnbauer 1995, Berkovits 1995, Beenakker et al 1995, Brouwer and Beenakker 
1995(a,b), Bruun et al 1995, Claughton et al 1995(b), Argman and Zee, 1996, Fraum et al 1996, Melsen et al 1996, 
Slevin et al 1996] 



VI. REENTRANT AND LONG-RANGE PROXIMITY EFFECTS. 



A. Reentrant behaviour of the conductance. 



As well as the ZBA discussed in the previous section, a number of recent experiments have revealed a non-monotonic 
behaviour of the voltage and temperature dependence of the conductance. In this section, we use the quasi-classical 
approach to highlight the origin of this phenomenon. For simplicity, we consider the case of a normal diffusive wire 
located in the region < x < L between a normal and a superconducting reservoir. As explained in section IV(C) 
the equation for the Green's function F R , which is conveniently written 

Dd 2 x u(x) + 2ie sinh(u(a;)) = 0, (145) 
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must be solved with boundary conditions u(0) = and u(L) = iir/2. The voltage dependent conductance is then 
obtained from 



G(e) 
Go 



— I dx 



1 



L J cosh 2 (Reu(x , e)) 



(146) 



where G is the conductance in the normal state. Using the identity cosh 2 w = 1 — tanh 2 u>, equation (146) 
becomes 



G(e) 
Go 



1 



= < 1 — — / dxtanh (Reu(x,e)) 



When Reu{x,e) is a small quantity, one can expand the denominator of equation (147) to yield 



5G(e 
G u 



1 f L 

— I dxtanh 2 (Reu(x,e)), 
L Jo 



(147) 



(148) 



which allows us to obtain approximate solutions in both the low- and high-energy limits. In the limit e <C ct where 
€t = D I L 2 , one can treat the term linear in the energy as a perturbation and write u(x, e) = u (x) + 5u(x, e). When 
inserted into equation (145), this yields 



u(x, e) 



7r x 8 e / . ,irx . x 

1- — — sin — ) 

2 L it 2 e T V y 2L J L 



(149) 



In the opposite limit of large energies, e> e-r, one notices that among the solutions of the original non- linear equation, 
there exists a subset satisfying the following equation 



d x u — 2k sinh(— ). 



(150) 



Solutions of equation (15C) are also solutions of the original equation (145) f| The solution of equation ( |15C| ) reads 

(151) 



tanh(^) = tanh(^M) e fc (^) 



where k = yj —2ie/ D. At the norm al reservoir, the value of the solution cannot be imposed, though in the limit 
of large energies, the r.h.s. of (151) is exponentially small and the solution satisfies the condition that the pairing 
function must vanish. We also note that in both the low- and high-energy limits, the real part of u(x, e) is a small 
quantity and therefore the assumption made at the beginning is justified. With these two solutions we may finally 
write 



G 



A 



(152) 



6G{£) =BM,e»e T , (153) 



Go V e 

with A and B numerical constants. 

The above analysis illustrates the key ingredients required to obtain a reentrant effect, which was noted prior to 
the experiments by [Artemenko, Volkov and Zaitsev 1979]. The effect originates from the presence of a distribution 
function which is spatially not in equilibrium and consequently will be most relevant in situations where the voltage 
drop is distributed along the wire, rather than dropping predominantly across the tunnel junction. 



5 Although the converse is not generally true, because there may be solutions of the original equation, which are not solutions 
of the above equation. This in general will depend on the boundary conditions, as is evident from the fact that the original 
equation is second order while the derived above equation is first order. 
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The fact that the conductance variation appears to be a small quantity both in the low energy and in the high 
energy limit, controlled respectively by the small parameters (e/er) 2 and y / er/e, is not accidental. It has been 
noticed that the conductance correction resembles the fluctuation correction in the theory of paraconductivity, the 
so-called Maki-Thompson term [Volkov 1994(a)]. In the high energy limit, the proximity effect induces Cooper pairs 
in the normal metal, which act like a superconducting fluctuation above the critical temperature. At low energy, the 
induced pair wave function is not small, so at first sight the correction should be sizeable. However the "longitudinal" 
component of the induced pair wave function gives a contribution which is equal to the normal state conductance and 
is not temperature dependent, while the "transverse" component, which is small at low energies, gives the additional 
temperature dependent contribution. On general grounds, one expects that the low energy e 2 dependence should be 
universal, while the high energy behaviour may differ from the inverse square root law, depending on the quality of 
the interface and on the effective dimensionality of the mesoscopic wire. 

B. Long-range proximity effects. 

We now discuss briefly the origin of the long-range proximity effects, observed in the experiments of [Petrashov et 
al. 1993,1995, Dimoulas et al. 1995, Courtois et al. 1996], where the conductance of a mesoscopic loop between two 
superconducting electrodes was measured. In such structures, the total current is the sum of two terms. The first is a 
Josephson current between the two superconducting electrodes and decays exponentially as L/Lt in agreement with 
known theory of superconducting weak links [Likharev 1979]. The second was shown to be due to a phase-periodic 
conductance with amplitude decaying as the inverse of the temperature. A detailed theory of the effect has been 
provided by [ Volkov, Allsopp and Lambert 1996 and Volkov and Takayanagi 1996(a)] using both quasi-classical and 
multiple scattering methods. The different behaviour of the two contributions to the current occurs because the 
Josephson current has a thermodynamic origin and is obtained by integrating over a range of energies of the order 
of fcfiT, while the second contribution giving rise to the reentrant behaviour is kinetic in nature. In the structure 
considered in the experiments, a non- vanishing value of the pair wave function F{x, e) is induced in the normal region. 
Away from the superconductor the pair wave function decays exponentially as 

F(x,e) w exp(-kx) (154) 

where k = \J e/D. The contribution to the Josephson critical current comes from an integral over the energy 
involving the product of anomalous Green's function F R F R and F A F A . Since this product contains Green's functions 
with poles on the same side of the real axis, one can deform the contour of integration and transform the integral 
over the energy to a sum over Matsubara frequencies e n = irT{2n + 1). As a result, the main contribution comes 
from energies of the order of the temperature. At a distance L, the Josephson current decays exponentially as L/Lt 
where Lt = tJD/T. By introducing the Thouless energy ct = D / L 2 , the condition for a vanishing Josephson current 
becomes er <C T. In addition to the above, the proximity effect gives rise to an additional contribution to the normal 
state conductance, which involves products of the type F R F A . As a consequence, the integral over the energy cannot 
be transformed to a sum over Matsubara frequencies. (This expresses the fact that the main contribution to the 
integral comes from energies of the order of the Thouless energy.) As a consequence, if the condition <C T is 
satisfied, then the pair wave function spreads over the entire sample. 

VII. ANDREEV INTERFEROMETERS 

When a quasi-particle Andreev reflects from a normal-superconducting interface, the phase of the outgoing ex- 
citation is shifted by the phase of the superconducting order parameter. Consequently if a phase-coherent normal 
conductor is in contact with two superconductors with order parameter phases fa and fa, transport properties will 
be oscillatory functions of the phase difference = fa — fa. Following a number of theoretical proposals [Spivak 
and Khemel'nitskii 1982, Altshuler and Spivak 1987, Nakano and Takayanagi 1991, Takagi 1992, Lambert 1993, Hui 
and Lambert 1993(a), Hekking and Nazarov 1993], several experimental realizations of Andreev interferometers have 
been reported [de Vegvar et al. 1994, Pothier et al. 1994, van Wees et al 1994, Dimoulas et al. 1995, Petrashov et 
al. 1995, Courtois et al 1996]. In addition to interferometers formed when a normal metal makes contact with two 
superconductors, the conductance of a single exended N-S structure is predicted to be an oscillatory function of the 
phase gradient across the interface [Cook et al 1995], although no experiments on such phase-gradiometers have been 
reported to-date. 
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These papers address a number of issues including the fundamental periodicity of the conductance oscillation, the 
nature of the zero-phase extremum, the magnitude of the oscillation amplitude and the role of disorder, geometry 
and dimensionality. In the main, experiments have probed either diffusive or ballistic structures, which in-turn can 
be divided into metallic samples with a conductance G much greater than e 2 /h and semiconducting structures with 
G less than or of order e 2 /h. For the former, the ensemble averaged conductance < G > is the relevant quantity 
and mesoscopic fluctuations are unimportant. These have confirmed the 2ir periodicity for < G > predicted by 
[Lambert 1993, Hui and Lambert 1993] and by all subsequent theories, [Hekking and Nazarov 1993, 1994, Zaitsev 
1994, Takane 1994, Allsopp et al 1996, Kadigrobov et al 1996, Nazarov and Stoof 1996, Volkov and Zaitsev 1996, 
Zagoskin et al 1996, Leadbeater and Lambert 1997]. For the latter, mesoscopic fluctuations can dominate and the 
non-averaged conductance G is of interest, which depends in detail on the geometry and impurity realisation of a 
particular sample. These confirm the 2ir perodicity for the non-averaged conductance predicted initially for diffusive 
systems by [Altshuler and Spivak 1987] and for one-dimensional clean systems by [Nakano and Takayanagi 1991 and 
Takagi 1992]. 

In this section we briefly review the various theoretical approaches and comment on their applicability to the 
experiments. 



A. Quasi-classical theory: the dirty limit. 

We start by using quasi-classical theory to compute the ensemble averaged conductance when the mesoscopic normal 
conductor is in the diffusive transport regime. As noted in previous sections, the influence of superconductivity on 
the transport properties of a mesoscopic conductor may change in a qualitative way, depending on the quality of the 
normal-superconductor interfaces. Consider first the zero temperature conductance of an interferometer comprising 
a tunnel junction (labelled as 1) connected by diffusive 1-D wires to a fork. Each of the two arms of the fork is a 
diffusive wire, connected via tunnel junctions (labelled as 2) to infinitely long superconductors. The conductance of the 
diffusive wires is assumed to be much greater than that of the tunnel junctions. In this structure the superconductors 
play the role of elctrode and the main voltage drop occurs at the tunnel barriers. The conductance can be analysed 
using the circuit theory outlined in section IV [Nazarov 1994, Zaitsev 1994], which yields 

G — ^G 2 G 2 cos 2 (0/2) 

{G 2 +4G 2 cos 2 (</>/2)} 3/2 

where Gi is the conductance of the tunnel junction (1), G2 is the conductance of the tunnel junction (2) and <f> is 
the phase difference between the two superconductors. In the limit G\ » G 2 , this simplifies to 

G = 4^1 cos 2 (0/2), (156) 



whereas if G 2 ^> G±, 



1 g| i_ 

2 G 2 I cos 2 (0/2) I' 



G - o TT „ c 2/ J ,/oM- ( 157 ) 



One can see that for G\ = G2 there is a zero-phase minimum, as there is for G\ <C G 2 . However when G\ 3> G 2 , 
there arises a zero-phase maximum. In all cases the conductance vanishes when the phase difference between the 
superconductors is n. Furthermore, the ratio Gi/G 2 controls the amplitude of the conductance oscillations which are 
greater when Gi w G 2 . 

The results of quasi-classical theory can of course be compared with the results of numerical multiple scattering 
theory If the parameters used in the numerical simulations are chosen in such a way to satisfy the assumptions 
underlying the quasi-classical treatment, then a detailed analysis shows that this is indeed the case [Claughton, 
Raimondi and Lambert 1996], although certain features, such as the vanishing of the conductance at 4> — ir, appear to 
be an artifact of te one-dimensional nature of the analysis leading to equation (155). It is perhaps worth emphasizing 
that the exact numerical techniques used in [Claughton, Raimondi and Lambert 1996] are not confined to the diffusive 
regime and an interesting result of these simulations is that when the resistance of a given structure is dominated 
by tunnel barriers, the diffusive nature of the wires is not relevant and similar results are obtained for both ballistic 
or diffusive wires. The geometry is clearly important, but when dissipation occurs along the wire, the almost one- 
dimensional approximation used in many quasi-classical calculations works well. 
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Another structure illustrating the role of finite temperature and voltage is suggested by the experiment by Pe- 
trashov et al. [1995], who measured the conductance of a mesoscopic wire, which makes contact with superconductors 
at points between the normal electrodes. At zero temperature, this structure is predicted to have only small am- 
plitude of osicllation arising from mesoscopic flcutuations, in contrast with the large amplitude oscillations observed 
experimentally. In the experiment however, a finite voltage drop is distributed along the wire and the temperature is 
non-zero. Consequently the phase periodic conductance is dominated by an effect similar to the reentrant behaviour 
discussed in the previous section. The analysis of this structure has been carried out in detail by [Nazarov and Stoof 
1996, Volkov, Allsopp and Lambert 1996, Stoof and Nazarov 1996], using both quasi-classical and numerical multiple 
scattering methods. 



B. Description of the ballistic limit: a simple two-channel model. 

Having discussed diffusive interferometers, we now examine the clean limit, which was first described in one dimen- 
sion by [Nakano and Takayanagi 1991 and Takagi 1992]. An analytic theory of the ballistic limit in higher dimensions 
has also been developed [Allsopp et al 1996], based on a multiple scattering description of a clean N-S interface. In the 
absence of disorder, when there is no phase difference between the two superconductors, translational invariance in the 
direction parallel to the interface (and transverse to the current flow) allows one to reduce the two-dimensional system 
to the sum over many independent one-dimensional channels. When a phase difference between the superconductors 
is imposed, inter-channel coupling is introduced, but as shown below, this coupling typically involves only pairs of 
channels in such a way that an accurate description is obtained by summing over independent pairs of coupled 
channels. This considerably simplifies evaluation of the boundary conductance 

N 

G B (0) = 2R a = 2Tr r a r\ = 2 £ {R a ) t] , (158) 

where (R a )ij — |(^a)ij| 2 is the Andreev reflection probability from channel j to channel i. Indeed the Andreev 
reflection coefficient is of the form R a = i?diag + -Roff-diag where i?diag = ^2iLi{Ra)u and -R ff-diag is the remaining 
contribution from inter-channel scattering, .Roff-diag = YliLj=i(Ra)ij • ^ channels only couple in pairs, both the off- 
diagonal scattering and the diagonal scattering will scale as the number of channels. Consider now a normal barrier 
in the N-region of a N-S interface. Particles (holes) impinging on the normal scatterer are described by a scattering 
matrix s pp , (shh), and those arriving at the N-S interface by a reflection matrix p, where 



o I PP pp I „ 

\ l pp r pp J \ Php Phh 

The elements of s and p are themselves matrices describing scattering between open channels of the external leads. 
For an ideal interface, where Andreev's approximation is valid, p pp and phh can be neglected and as a consequence, 
Php and p p h are unitary and one obtains a generalization of a formula due to Beennakker [1992] r a — t' hh phpM pp t ppi 
with M pp = 1 — r' pp p p hr hh php. In contrast with the analysis of Bennakkcr [1992], where php is proportional to the 
unit matrix, the interference effect of interest here is contained in the fact that ph p induces off-diagonal scattering. 



tpp \ n — I Ppp Pp h 



Substituting r a into equation (158) and taking advantage of particle-hole symmetry at E = 0, yields 

G = 2Tr(TQ- 1 T(Q' 1 )- 1 ) (159) 

where Q = p ph + {r')ppPph{r')pp, with T = t pp t pp the transmission matrix of the normal scattering region. This 
multiple scattering formula for the boundary conductance is valid in the presence of an arbitrary number of channels 
and in any dimension. 



Equation (159) is very general and makes no assumption about the nature of matrices p p h and s pp . In a two-channel 



model, p p h is chosen to be an arbitrary two dimensional unitary matrix and since in the absence of disorder, t pp and 



Tpp are diagonal and interchannel coupling arises from p p h only. Substituting these matrices into equation (159), yields 
an expression for r a involving a single phase whose value is a linear combination of phase shifts due to normal 
reflection at the barrier, Andreev reflection at the N-S interface and the phase accumulated by an excitation travelling 
from the barrier to the interface. After averaging over the rapidly varying phase 9, the diagonal and off-diagonal 
terms are determined and found to vary periodically with with period 2ir as predicted in one dimension by [Nakano 
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and Takayanagi 1991] A key new feature of equation ( |159| ) is that the in the absence of the normal barrier, unitarity 
ensures the the total Andreev reflection obtained by summing these two contribution is independent of <fi and therefore 
the amplitude of conductance oscillation vanishes. Thus to obtain a maximum amplitude of oscillation, the tunnel 
barrier must be tuned to a non-zero value. 

An interesting mechanism capable of producing large conductance oscillations which scale with the number of 
channels has been proposed by [Kadigrobov et al 1995], based on the simultaneous resonance of Andreev levels 
corresponding to different channels. To illustrate this, consider for example the clean S-N-S structure analyzed by 
[Bardeen and Johnson 1972 and Kulik 1970]. For a long-junction, the Andreev levels have energies 

£ n = ^HV±(2n + l)7r], (160) 

where vf,z is the component of the Fermi velocity perpendicular to the N-S interfaces and A</> the phase difference 
between the superconductors. In general these energies do not coincide, but when Acj) = (2m + l)n all levels pass 
through the Fermi energy E — 0, yielding large-scale oscillations in transport properties. In the calculation of 
[Kadigrobov et al 1995], normal scattering arises from the beam splitters which form an integral part of the structure 
considered and therefore no addition tunnel barriers are needed to break the unitarity condition of [Allsopp et al 
1996]. 



VIII. CONDUCTANCE SUPPRESSION BY SUPERCONDUCTIVITY. 



When a superconducting island is added to a normal host, one naively expects that the electrical conductance 
of the composite material will increase. In contrast detailed calculations based on the multiple scattering formula 
[Hui and Lambert 1993(b), Claughton et al 1995] predicted that the change 8G in the two probe conductance of 
a mesoscopic sample, due to the onset of superconductivity can have arbitrary sign. Indeed for a sample with a 
high enough conductance, a simple theorem [Hui and Lambert 1993] states that 5G is guaranteed to be negative 
and furthermore the magnitude of SG is predicted to scale with G. This "anomalous proximity effect" occurs in the 
absence of tunnel barriers and at first sight appears to conflict with the quasi-classical description presented in section 
IV, which predicts that at zero temperature, the conductance of equation (27), which takes the form G = (4e 2 /h)R a 
is unchanged by the onset of superconductivity. In this section we briefly discuss the origin of this effect and how it 
can be reconciled with quasi-classical theory. 

For convenience we separately discuss the clean, diffusive and almost Anderson localised limits. 



A. The clean limit. 



The Landauer formula (2) predicts that the conductance of a normal ballistic wire is G = (2e 2 /h)N, where N is the 
number of open scattering channels at the Fermi energy Ep. Consequently in the absence of spin-splitting, if Ep is 
varied, G exhibits a series of steps of height 2e 2 /h corresponding to the opening or closing of scattering channels. The 
effect of inducing a uniform superconducting order parameter of magnitude Ao in such a wire is shown in figure 9 of 
[Claughton et al 1995a]. In the vicinity of normal-state steps, the conductance is predicted to be extremely sensitive 
to the onset of superconductivity. With increasing Ao, the steps are destroyed and G decreases. This behaviour is 
attributed to the breakdown of Andreev's approximation, which is most pronounced for scattering channels whose 
wavevector k x normal to the N-S interface tends to zero. Such channels occur precisely at a normal-state conductance 
step. 

For a clean N-S-N structure, provided the superconductor is much longer than the superconducting coherence 
length £ = kp 1 Ep / Aq, there is negligible quasi-particle transmission through the S region and therefore the total 
resistance reduces to the sum of two boundary resistances. Since the boundaries are identical, this yields for the 
total conductance defined in equation (|30|), G = R a and since the system consists of decoupled channels, R a can 
be obtained by solving the Bogoliubov - de Gennes equation at a one dimensional N-S interface. By insisting that 
scattered wavefunctions and their first derivatives be continuous at the boundary, one finds [Claughton et al 1995a] 

N 

R a = J2 R ^ ( 161 ) 

where 



29 



and iJ, n = h 2 (k^) 2 /2m is the longitudinal kinetic energy of a quasi-particle incident along channel n. Clearly those 
channels corresponding to low angle quasi-particles with \x n < A , possess a small Andreev reflection probability R n 
and since there is no transmission, the corresponding normal reflection probability (1 — R n ) approaches unity. Finally 
one obtains for conductance change 5G due to the switching on of a uniform order parameter in such a clean system, 



It is perhaps worth noting that although the suppression of Andreev reflection for low-angle quasi-particles is quite 
general, the destruction or otherwise of conductance steps in a quantum point contact, depends on the geometry 
of the contact. If quasi-particles are adiabatically accelerated before encountering the N-S interface, then Andreev 
reflection need not be suppressed and conductance steps can survive. The criterion for the survival of the nth step is 
clearly 



The results of the previous subsection reveal that 5G is sensitive to fine scale structure in the normal state conduc- 
tance Gm- For clean systems this structure takes the form of well-known conductance steps. For strongly disordered 
systems, it is known that Gjv can exhibit sharp resonances and therefore the behaviour of G with increasing Ao 
is sensitive to such features. This behaviour is illustrated in figure 10 of [Claughton et al 1995], which shows that 
negative changes in the conductance with increasing A are associated with resonances in the normal-state transport. 

In view of tunnelling theory, one might regard the suppression of conductance in an almost insulating N-S structure 
as unsurprising, in which case the increase in the conductance due to the onset of superconductivity would be 
regarded as anomalous. A detailed description of such resonant transport in N-S and N-S-N structures and in 
resonant interferometers is provided by [Claughton and Lambert 1995b], which follows an earlier zero- voltage theory 
[Bccnakkcr 1992] and a one-dimensional description of a delta- like potential well [Khlus et al 1994]. 



The above results show that the conductance of a clean or very dirty conductor can either increase of decrease 
when superconductivity is induced, depending on the microscopic impurity configuration and on the geometry. This 
qualitative behaviour has been observed in experiments by Petrashov and Antonov [1991] and Petrashov et al [1993(b)] , 
which also exhibit the weak magnetic field behaviour shown in fig 12 of [Claughton et al 1995a]. Nevertheless a 
quantitative comparison between these calculations and experiments is not possible, because of the differing geometries 
and multiprobe measurements. 

A first attempt at a quantitative comparison has been made recently by Wilhclm, Zaikin and Courtois [1997] 
who noted that although the one-dimensional approximation used in sec. V and VI to study transport properties 
of N-S structures successfully explains a range of experiments, the experiment by [Petrashov et al 1991] requires a 
multiprobe description. The key observation is that the kinetic properties of a two-dimensional metallic film, in contact 
with a superconductor, differ substantially from those of a quasi one-dimensional wire, because of an inhomogenous 
distribution of the currents in the sample. Although a detailed explanation of the effect can be obtained by solving 
the equations for the quasi-classical Green's function in a two-dimensional structure, it is possible to gain considerable 
insight by adopting an effective circuit model, which captures the main features. For example consider four wires 
arranged in a square. The top two corners of the square (labelled A and B) are the voltage probes, while the lower 
two corners (labelled C and D) are the current probes. Let the left-most vertical wire between corners A and C be 
placed in contact with a superconductor and labelled 1 , the right-most vertical wire between corners B and D labelled 
2 and the upper and lower wires labelled 3 and 4 respectively. Then the four-probe conductance reads 




(163) 



A /fi n « 1. 



(164) 



B. The localised limit and resonant transport. 



C. The diffusive limit. 



G four — G3,Gi(G 1 + G 2 + G 3 + G 4 ), 



(165) 
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where Gi is the conductance of wire i. The analysis of sections IV, V and VI shows that the conductance of a wire 
in contact with a superconductor is changed in several ways. 1) The local conductivity is increased by the presence 
of a superconductor, but this increase decays over a distance from the superconductor of the order of the coherence 
length of the normal metal D/max{T, E}. 2) The conductance has a non monotonic behaviour as function of the 
energy and temperature. To a first approximation one may consider that when T 3> D/ L 2 , only the wire attached 
to the superconductor will be appreciably affected by the presence of the superconductor. In such a situation the 
conductances G2, G3 and G4 are the same as in the normal state, while G\ is increased. In this case, according to 
equation (165), the four-probe conductance is decreased by superconductivity. On the other hand, at T -C D/L 2 , all 
conductances increase and so does the four-probes conductance. As [Wilhelm ct al 1997] have shown this leads to a 
much richer structure in the temperature and voltage dependence of a two-dimensional film, compared with that of 
a one-dimensional wire and to qualitative agreement with experiment. 

A key prediction of [Wilhelm et al 1997] is that at zero bias and zero temperature this effect should vanish and 
the conductance should be unchanged by the onset of superconductivity. In contrast [Seviour et al 1997] identify an 
additional large-scale effect due to weak localisation, which is not contained in standard quasi-classical theory and 
which will survive at zero-temperature and zero-bias. 



IX. CONCLUSION. 



We conclude this review by indicating a few lines of future research. The results of the first experiments in 
hybrid structures have been successfully explained in terms of almost one-dimensional models. This has been possible 
because the physical phenomena responsible for the transport properties of a given structure mainly occurred at 
the N-S interface, so that the actual geometrical arrangement of current and voltage probes was not important. 
This is clearly the case in the presence of low transmitting interfaces, which tend to dominate the overall transport 
properties. The improvement of the quality of N-S interfaces betweens separate parts of the hybrid structure and 
the contact leads makes necessary to take into account the detailed geometrical arrangement of a given measurement 
setup [Allsopp, Hui, Lambert and Robinson 1994]. Under this respect, the last twelve months have witnessed a 
considerable experimental effort at paying particular attention to the geometry of a given measurement. As we have 
seen in the previous section taking into account the actual disposition of the current and voltage probes has been 
crucial to describe the effect of the conductance suppression, even though the structure itself has a one-dimensional 
geometry [Wilheml, Zaikin and Courtois 1997]. The recent experiment by [Hartog et al 1996] has beautifully shown 
the emergence of new physical effects when the geometry of the experimental setup is considered carefully. 

A second important semplifying assumption in examining the experiments in hybrid structures consisted in ne- 
glecting the effect of electron-electron Coulomb interaction. This is a reasonable assumption in many metallic and 
semiconducting wires used in the experiments. In these cases, a mean-field treatment of electron-electron interaction 
is sufficient and can be easily incorporated in the effective electronic parameters. In the last decade there has been 
however a continuous progress in downsizing fabrication of electronic systems, so that it is now possible to measure 
transport through almost zero-dimensional structures, or quantum dots. In these systems, the effect of Coulomb 
interaction is no longer negligeable, and various new effects arise. Besides quantum dots, quantum wires also have 
been a subject of noticeable attention. It is well known that in one-dimensional systems electron-electron interaction 
drastically modifies the low energy properties and lead to the concept of Luttinger liquid. Quantum wires are believed 
to be an experimental realization of the Luttinger liquid state. As a consequence in the last few years a considerable 
theoretical literature has accumulated with the intent of studing phase coherent transport in interacting systems. It 
is almost a logical consequence to foresee in the forthcoming years the emergence of a new field of research aimed 
at combining the fields of mesoscopic superconductivity and strongly interacting systems. Indees already quite a few 
papers have started to appear [Fisher 1994, Fazio, Hekking and Odintsov 1995, Maslov, Stone, Goldbart and Loss 
1996, Fazio and Raimondi 1997, Takane and Koyama 1997] and the time is mature for new exciting developments. 
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FIG. 1. Various generic experimental arrangements used in measuring dc transport in superconducting hybrids. The grey 
area indicates a phase-coherent normal region, whereas superconductive parts are marked by a S. The widening open parts at 
the ends represent the reservoirs. 
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